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[57] ABSTRACT 

A digital phase-lock loop (DPLL) which generates a 
signal with a phase that approximates the phase of a 
received signal with a linear estimator. The effect of a 
complication associated with non-zero transport delays 
related to DPLL mechanization is then compensated by 
a predictor. The estimator provides recursive estimates 
of phase, frequency, and higher order derivatives, while 
the predictor compensates for transport lag inherent in 
the loop. 

3 Claims, 9 Drawing Sheets 
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Phase-Locked Loops”, IEEE Transactions on Commu- 
DIGITAL PHASE-LOCK LOOP HAVING AN nications Technology, Vol 16, No. 2, April 1968, at- 

ESITMATOR AND PREDICTOR OF ERROR tempts to optimize a concise performance measure, such 

as a weighted sum of integrated transient response, gain 
ORIGIN OF INVENTION 5 margin, and noise variance. For some cases, analytical 

The invention described herein was made in the per- expressions for the performance measure can be ob- 

formance of work under a NASA contract, and is sub- tained, and closed-form expressions for loop filters de- 
ject to the provisions of Public Law 96-517 (35 USC fined. Disadvantages are that the loop filter form (e.g., 

202) in which the Contractor has elected not to retain order of loop) is not defined apriori, but rather is a result 

title. 10 of the optimization procedure, and that the design pro- 


BACKGROUND OF THE INVENTION 

This invention relates to a digital phase-lock loop 
(DPLL), and more particularly to a DPLL in which the 
transfer function of the loop is comprised of an estima- 15 
tor and a predictor of error. * 

A phase locked loop (PLL) is a feedback system that 
tracks the phase of a received quasi-periodic signal. A 
typical PLL structure is shown in FIG. 1. The differ- 
ence between the received phase and a voltage con- 20 
trolled oscillator (V CO) phase is filtered and used to so 
adjust the VCO that the VCO phase tracks the received 
phase. 

Designers of PLL’s select a loop filter to meet perfor- 
mance requirements, usually specified in terms of band- 25 
width, gain margin, and dynamic errors. Design meth- 
ods are well documented in the literature. Gardner, F. 

M., Phaselock Techniques , John Wiley and Sons, (1979). 
Best, R. E., Phase Locked Loops , Theory. , Design and 
Applications , McGraw Hill, (1984). Viterbi, A. J., Princi - 30 
pies of Coherent Communication , McGraw Hill (1966). 

What is of interest here is a digital phase-lock loop 
(DPPL) characterized by constant loop update rates. 
For a discussion of DPLLs with variable loop update 
rate, see Gill, G. S., and Gupta, S. C., “First Order 35 
Discrete Phase Locked Loops with Applications to 
Demodulation of Angle Modulated Carrier”, IEEE 
Trans, on Communications, June 1972; Weinberg, A. 
and Liu, B., “Discrete Time Analysis of Nonuniform 
Sampling First-and-Second-Order Digital Phase 40 
Locked Loops”, IEEE Trans, on Communications, 
February 1974; Osborne, H. C., “Stability Analysis of 
Nth Power Digital Phase Locked Loop — Part I: First 
Order DPLL”, IEEE Trans, on Communications, Au- 
gust 1980; Osborne, H. C., “Stability Analysis of Nth 45 
Power Digital Phase Locked Loop — Part II: Second 
and Third Order DPLLs”, IEEE Trans, on Communi- 
cations, August 1980. DPPLs with fixed update rate are 
often designed either by analogy to continuous time 
domain PLLs, or based on an optimality criterion. 50 

Design by analogy to continuous domain PLLs, i.e., 
applying s-plane design rules by translation to the z- 
plane, suffers from two major disadvantages. The de- 
sign is based on iteratively placing open-loop poles and 
zeroes at “well chosen” locations until satisfactory per- 55 
formance is achieved. Such locations are selected based 
on design experience rather than on well established set 
of rules. Also, the design does not accound for the trans- 
port delay in the digital loop, caused by hardware de- 
lays and loop filter computations. This delay is either 60 
“tolerated,” i.e., the degradation due to the delay is 
analyzed and found to be acceptable or “compensated” 
by the equivalent of a lead-lag network. This design 
procedure is usually acceptable only when the loop 
update rate is very high compared to the loop band- 65 
width. 

Design based on an optimality criterion as discussed, 
for example, by Gupta, S. C., “On Optimum Digital 


cedure is complex. 

The estimator-predictor approach described in this 
application offers an alternate design for DPLLs. The 
DPLL design presented here is shown to be equivalent 
to an estimator followed by a predictor compensating 
for the transport lag. The estimator incorporates past 
measurements of phase into a state vector estimate, 
consisting of phase, frequency, and perhaps higher de- 
rivatives. Selection of an estimator can utilize the base 
resources available in linear optimal estimation theory. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, the basic 
purpose of a DPLL (to generate a signal with a phase 
that approximates the phase of a received signal) is 
accomplished with a linear estimator. The effect of a 
complication associated with non-zero transport delays 
related to DPLL mechanization is then compensated by 
a predictor. The estimator provides recursive estimates 
of phase, frequency, and higher order derivatives, while 
the predictor compensates for the transport lag inherent 
in the loop. The closed-loop DPLL transfer function 
has the form 

H(z)=C(z)D(z)z-" 

where H(z) is the closed-loop DPLL transfer function, 
C(z) is the estimator transfer function, D(z) is the pre- 
dictor function, and z-^is the DPLL transport lag in 
units of loop update time. A phase detector with gain A 
is used to implement this equation in a closed-door 
DPLL. Thus, the transfer function H(z) is given by the 
above equation only where the phase detector has a 
nominal gain value of A = 1. One realization of the esti- 
mator, predictor and delay is comprised of a digital 
filter having a transfer function S(z) and a numerical 
control oscillator having a transfer function Q(z), where 
S(z) satisfies. 



In a preferred embodiment, a detector and a numerical 
control oscillator are implemented with hardware in a 
direct connection to the digital filter implemented in a 
programmed digital computer in which the function 
Q(z) of a numerical control oscillator is implemented as 
part of the computer program. This decomposition of 
the DPLL closed-loop transfer function results in a 
straightforward loop filter design, enabling use of tech- 
niques from optimal and sub-optimal estimation theory. 

BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 illustrates a prior-art phase lock loop (PLL). 

FIG. 2 illustrates a functional block diagram of a 
digital phase-lock loop organized in accordance with 
the present invention. 
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FIG. 3 is a block diagram which may be readily im- 
plemented from the functional block diagram of FIG. 2. 

FIG. 4 illustrates a hybrid software and hardware 
implementation of the DPLL shown in FIG. 3. 

FIG. 5 is a graph of dynamic phase error versus esti- 5 
mator parameters. 

FIG. 6 is a graph of loop bandwidth versus estimator 
parameters. 

FIGS, la to Id illustrate examples of root loci for: 

(a) Second order estimator; 10 

(b) Third order estimator; 

(c) Second order estimator, with pole at z— —0.9; 
and 

(d) Third order estimator, with pole at z = — 0.9. 

FIG. 8 is a graph of upper gain margin versus normal- 15 

ized loop bandwidth.’ 

FIG. 9 is a graph of lower gain margin versus normal- 
ized loop bandwidth 

FIG. 10 is a graph of phase error coefficients versus 
normalized loop bandwidth. 20 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

The basicjpurpose of a DPLL is to generate a signal 
with phase 0(z) that approximates 0(z), the phase of a 2 
received signal. From the viewpoint of linear estimation 
theory, this is an estimation problem, and the solution is 
found by the following procedure. First a linear state 
model for 0(z) is defined, with a corresponding mea- 3Q 
surement model. Then, statistical models for state (or 
process) noise and measurement noise are developed. 
Finally, some type of optimal estimator is selected such 
as a recursive least squares estimator. 

There is a slight complication since a non-zero trans- 35 
port delay is associated with DPLL mechanization. The 
effect of this delay is compensated by a predictor that 
extrapolates the state estimate by an appropriate time 
interval. Thus, the closed-loop DPLL transfer function 
has the form: an 


S(z) = 


c&m 

Qiz) 


z-X 


( 2 ) 


FIG. 4 illustrates a preferred DPLL implementation. 
The phase detector 13 and “hardware” NCO 17a are 
implemented with hardware, while the summer 15, 
digital filter 16, and a “software” NCO lib are imple- 
mented with a programmed computer. Both “hard- 
ware” NCO and “software” NCO have identical math- 
ematical representations, as the latter simulates the op- 
eration of the former. When comparing FIGS. 1 and 4, 
it is apparent that the function of a loop filter in a PLL 
of FIG. 1 is performed by a minor loop 20 in FIG. 4 
implemented in the computer-resident portion of the 
DPLL. 

The design of a DPLL conforming to FIGS. 2, 3 and 
4 can be divided into five distinct steps. The first step is 
to select a model for the received phase process. This 
model can include the process dynamics and statistics 
on process noises and measurement noise. The second 
step is to model the “hardware” NCO 17a, transport lag 
12, and phase detector 13. As shown below, the trans- 
port lag may be incorporated as part of the transfer 
function Q(z) of the NCO lib. 

The next step is selection of an estimator 10, based on 
the phase model and the process noise and measurement 
noise covariance matrices. Possible realizations of the 
estimator are some variations of Kalman filters, or other 
least squares estimators. The fourth step is determina- 
tion of a predictor to compensate for transport lag. The 
last, but definitely not the least, step is to assure that the 
loop is stable. 

Two separate models are used for received phase. 
The first model has two state variables: phase and fre- 
quency, while the second model uses three-state vari- 
ables: phase, frequecy and frequency rate. These models 
are addressed below as second-and third-order models, 
respectively. The time-invariant phase models are ex- 
pressed as: 


H(z)=C(z)D(z)z- n 


( 1 ) 


X n — <Px n — I -j- _1 


(3) 


where H(z) is the closed-loop DPLL transfer function, 
C(z) is the estimator transfer function, D(z) is the pre- 
dictor transfer function, and z ~ N is the DPLL transport 
lag in units of loop update time. It is assumed here that 
the DPLL is updated at uniformly spaced time points, 
with update interval of T, hence the transport lag is an 
integer multiple of T. 

FIG. 2 shows a DPLL functional block diagram that 
addresses eq. (1). The three functional blocks 10, 11 and 
12 on the right are the elements of H(z), namely an 
estimator C(z), a predictor D(z), and an inherent trans- 
port lag z~ N . In addition, the figure incorporates a 
phase detector 13 with gain A represented by block 14 
and a summer 15. Note that when the phase detector 
gain has a nominal value of A= 1, the feedback is can- 
celled, and the closed loop transfer function is H(z) 
given by eq. (1). The closed loop transfer function, for 
general phase detector gain, is denoted P(z). 

The functional diagram shown in FIG. 2 is readily 
implemented in accordance with a block diagram pres- 
ented in FIG. 3. The estimator-predictor-delay combi- 
nation is replaced by actual components: a digital filter 
16 having a transfer function S(z), and a numerically 
controlled oscillator (NCO) 17 having a transfer func- 
tion Q(z), where S(z) satisfies: 


yn—H x n+ u n (4) 

where x n is the state vector, y n is the measurement, <P is 
the state transition matrix, u w is the measurement noise, 
and Vfl_i is the state, or process noise. Expressions for 
<I>, H, and the measurement covariance, R, are pres- 
ented in Table 1. Since only phase is measured, R is a 
scalar. 


TABLE 1 



Matrices for time invariant signal model 


Second order 

Third order 

<t> 

(«) 

( I T T 2 /2 ^ 
01 T 
\oo i ; 

H 

(1,0) 

0,o,o) 

R 

(<r„ 2 ) 

{CTu 2 ) 


Model of Phase Detector and NCO 

The phase detector inputs are actually time continu- 
ous signals. The phase detector measures the average 
phase difference between the phase of the input process 
and the NCO phase. This is equivalent to the difference 
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between the average input phase and the average NCO 
phase. Thus, 0( z) represents the input phase, averaged 
over the measurement interval, and 0(z) represents the 
NCO phase, averaged over the same interval. 

Hurd and Aguirre, and Simon and Mileant, cited 5 
above, discuss a general form for a DPLL NCO. The 
NCO is updated gT seconds after the phase measure- 
ment is made, where T is the loop update interval and 
0<— g«— 1. The case g=0 corresponds to NCO update 
immediately following the measurement, while g= 1 is 10 
caused by a T seconds delay before the NCO update. 
Such a delay occurs in practical systems while the com- 
puter determines the next NCO input. The transfer 
function Q(z) is the transfer function between the NCO 
input, and the average NCO phase at the next measure- 15 
ment, including all delays. This transfer function is: 


Q, } J, aoz) 2 ^ + (1 4- 2z — l^z 4- g 2 
^ 2 z2(* - 1) 

(5) 

20 

In particular, for the cases g=0 and g= 
the simple forms: 

1, Q(z) has 

* = ° ow-f- 

(6) 25 

g = i m = f^-j~ 

(7) 


6 

solution for equations 3-4, derived in Statmen, cited 
above, for a specific form of the process noise covari- 
ance matrix is: 

x nn — + K{y n — i f «— i) (9) 

where R -1 is the measurement noise covariance matrix 
and M is the solution for the linear matrix equation: 

M = H t R~ x H + (10) 

a = e^/r 

The parameter r is denoted the filter time constant. 
Using these equations, a z-plane filter transfer function 
can be evaluated. For a filter with m inputs (measure- 
ments) and n outputs (states), the transfer function is a 
matrix C(z) with n rows and m columns. C(z) can be 
computed from: 

C(z)=(zI-4>+KH<P)- l Kz (11) 

Table 2 summarizes the elements of the transfer func- 
tion. Since only phase measurements are used, the trans- 
fer functions are column vectors. Note that all the trans- 
fer functions poles are at z = — a within the unit circle, 
hence assuring filter stability. 


TABLE 2 



Transfer functions from input phase to state parameter 

State parameter 

Second order 

Third order 

Phase 

(t - a 2 )z ( z - 2 T+"a* ) 

(1 — a)z((a 2 -1- a 4- l)z 2 — 3a(l 4- a)z 4- 3a 2 ) 
(z - a) 3 


(z - a) 2 


Frequency 

(l - a) 2 z(z - 1) 

(1 ~ a)z(z - M3 - 3a 2 )z + 5a 2 - 4a - 1) 


H 

'n' 

1 

4 , 

2T(z - a) 3 

Acceleration 

— 

O - a) 3 z<z - l) 2 
T 2 (z — a) 3 


Note that in equations (6) and (7), the NCO function 
can be separated into two parts: integer number of de- 
lays, and a kemal, Qi(z). Qi(z) is characterized by nu- 
merator and denominator of the same degree. The cases 
of g=0 and g— 1 are then examples of N= 1 and N=2 
computational delays, respectively. The general expres- 
sions for a delay of N update times is: 

®=f f^f (8) 

In particular, cases with N=3, 4 . * . corresponds to 55 
pipelined feedback computations. 

Selection of an Estimator 

For this example, a least squares, fading memory 
estimator, with fixed gains, is selected. Such an estima- 60 
tor, described by Statman, J., “Simplified Solution for a 
Class of Fading Memory Filters,” submitted to IEEE 
Transactions on Aerospace and Electronic Systems, 
computes a state estimate based on a set of measure- 
ments, applying an exponentially decaying weight to 65 
past data. This “aging function” effectively discards 
measurements that are older than three or four estima- 
tor time constants. The steady state fading memory 


Selection of a Predictor 

The predictor generates a phase predict by extrapo- 
lating from the current phase estimate using either an 
estimate of frequency for a second-order model, or 
50 estimates of frequency and frequency rate for a third- 
order model. The extrapolation time, NT, is the number 
of integer delays in the loop multiplied by the update 
interval. Predictor equations are given in Table 3. 


TABLE 3 



Predictor formulas 
Second order 

Third order 

D(z) 

(1,NT) . 

1,NT, -QTO 2 ) 


( 


Loop Transfer Functions 

Recall that the transfer function H(z) equals 
D(z)C(z)z~ iV . Thus, in the example, H(z) depends only 
on the order of the phase model, the number of delays, 
and the estimator decay factor, a. The open loop, 
closed loop, and error transfer functions (G(z), P(z) and 
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E(z), respectively) for the DPLL, are easily expressed 
in terms of the function Ei(z)= 1 — H(z): 


, e(z) AH(z) A ( l - £ lW) 

- <t>(z) - 1 - Hiz) - E X {z) 

(12) 

5 

J(z) Giz) AH(z ) 

B(z) ~ 1 + G(z) ~ 1 4- (A - 1 )H(z) ~ 

03) 


A{\ - E X (z)) 
l + (A - 1)(1 - £i(z)) 10 


8 

a zero at z= — 1, thus, S(z) has a pole at that point. This 
means that S(z), were it not a part of the feedback loop, 
would be marginally stable. H(z) is stable since the pole 
at z = — 1, contributed by S(z), is cancelled by the zero 
5 of Q(z) at the same location. To maintain numerical 
stability, this pole of S(z) may be moved slightly into the 
unit circle. The next section demonstrates that this shift 
in pole location has negligible effect on loop perfor- 
mance. 

Performance Analysis 


m = %§- = i-m = 


E\(z) 


1 + (4 - 1)0 - E X {z)) 


where <I>(z) is the z-transform of the phase out of the 
phase detector, not including phase detector gain, i.e., 
<!>(z)=0(z)— (0(z). Expressions for Ei(z) are presented 


(14) In this section performance of the loops is analyzed. 
First, loop type and error due to dynamics of the re- 
ceived phase are evaluated, then loop bandwidths are 
15 determined. And finally, gain margins are computed. 
The discussions of dynamic error and loop bandwidth 
are for the nominal operating point, A=l. 


in Table 4. 

TABLE 4 


Formulas for E;(z) 

N Second order 

Third order 

1 - n 2 

<2 ~ D 3 

(z-a) 2 

(z - a) 3 

2 (z- 1)2 z + 2(1 - a) 

(z - l) 3 z + 3(1 - a) 

(z - a) 2 z 

(z — a) 3 z 

3 (z — l) 2 z 2 4- 2(1 — a) z 4- (a 2 

- 4a 4- 3) (z - l) 3 z 2 + 3(1 - a)z 4- (3a 2 - 9a 4- 6) 

(z — a) 2 z 2 

(z - a) 3 z 2 


A special case of interest occurs when the phase de- 
tector gain, A, has its nominal value A= 1. Then: 

G(z) = P(z) = H(z), E(z) = Ei(z) (15) 


Realizability and Stability 

The DPLL structure of FIG. 4 has three blocks that 
need to be realized: a “hardware” NCO, a “software” 
NCO and a digital filter S(z). The “hardware” NCO is 
realizable by definition. The “software” NCO is a mere 
simulation of the “hardware” NCO function and is thus 
realizable. This leaves the digital filter S(z). A necessary 
and sufficient condition for the realizability of this filter 
is that the order of the numerator is less or equal to that 
of the denominator. Oppenheim, A. V. and Schafer, R. 
W. “Digital Signal Processing”, Prentice-Hall, N.J. 
1975. Recall, from equations (2) and (8) that 
S(z)=D(z)C(z)/Qi(z). Since the above condition is 
satisfied for each component of C(z), it is also true for 
D(z)C(z) (the predictor merely computes a fixed weight 
sum of components of C(z)). The above condition is also 
satisfied for 1/Qi(z) (eq. (8)). Thus, S(z) is realizable. 

The following restricts the stability discussion to the 
nominal operating point, A=l, with analysis of the 
range of stable operating conditions deferred to the next 
section. When A=l, the closed loop function is: 

I\z)=H(z)— 1 — i?i(z) (16) 

Equation (16), combined with Table 4, shows that 
P(z) has poles at z=a and at z=0. Since a is real and 
satisfies 0<a< 1, all the poles of P(z) are inside the unit 
circle and the closed loop transfer function is stable. 
Gain margin is computed in the next section. 

Even though P(z) is stable, special care must be exer- 
cised in the implementation of S(z). Note that Qi(z) has 


(1) Loop Type and Dynamic Tracking Errors 

Loop type is the number of integrators in the open 
loop transfer function, i.e., the number of poles at z= 1. 
From equation (12), observe that the poles of the open 
loop transfer function are the same as the zeros of Ei(z). 
Table 4 presents Ei(z) for the second order and third 
order phase models, showing that Ei(z) has as many 
zeros at z=l as the order of the estimator, i.e., second 
and third order phase models correspond to type II and 
type III loops, respectively. Note that while loop type is 
independent of the number of computational delays, the 
order of the transfer functions strongly depends on the 
number of delays. Specifically, the orders of the numer- 
ator and denominator of Ei(z) are equal to the loop type 
plus N— 1, where N is the number of delays. 

The steady state phase error for constant dynamics is 
computed using the final value theorem: 

(17) 

E a = Hm 2 1 ~~ * E(z)6(z) 

2-1 2 

A type II loop has zero phase error for phase ramp, 
constant (non zero) phase error for phase acceleration, 
and infinite phase error for phase jerk, or higher dynam- 
ics. In a similar way, a type III loop has zero, constant, 
and infinite phase errors for phase acceleration, jerk, 
and higher dynamics, respectively. Thus the dynamics 
of interest are phase acceleration for a type II loop, and 
phase jerk for type III loop. Table 5 lists the steady state 
phase errors for the types II and III loops, at the operat- 
ing point, A=l. FIG. 5 shows these phase errors as 
function of the estimator time constant. The inputs are 
unit phase acceleration, e.g., rad/s 2 or m/s 2 , for the type 
II loop and unit phase jerk for the type III loop, and the 
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errors are in consistent units of rad or m, respectively. 
Note that as the estimator time constant increases, the 
phase error increases. In fact, for T < < r and N= 1, the 
normalized phase errors are r 2 and r 3 for the type II and 
type III loops, respectively. The phase error also in- 5 
creases slightly with the increase in the number of de- 
lays. This is caused by the increased prediction error. 

TABLE 5 

Steady state phase errors due to dynamics (phase errors) 

Second Order Estimator Third Order Estimator 

n (w> = T‘ 2 ) ( ew = T‘ 3 ) 


1 -p 2 

(1 - a) 2 

2 -p2 


(3 - 2a) 


(1 - a) 3 
T 3 


(1 - a) 3 (l-a) 3 

3 t 2 t T 3 

* — (a 2 - 6a + 6) k 


(1 - a) 2 


0 - a) 3 


(4 - 3a) 


(3a 2 - 12a + 10) 


15 


20 


In using the estimator-predictor formulation, it is 
important to use a good predictor. For example, if one 
uses the current phase as the predict, i.e., 

D(z)={ 1,0) or Z)(z)=(l,0,0) (18) 


then both the second order and third order estimators 
result in type I loops. 

(2) Bandwidth Results 

The one-sided bandwidth of the loop is defined as: 


o i j. L_ 

IT [P(l )2 277 j 


unit Z 


09) 


35 


where P(z) is the closed loop transfer function, and the 
integral is on the unit circle. At the desired operating 
point, A= 1, P(z)=H(z), and the integral can be written 
as: 


1 ( 20 ) 
B L T=k f \H(^)\2dx 
0 

The quantity B^T is also called the normalized band- 
width. Equation (20) uses the fact that in the example, at 
z=l, H(z) is unity, as can be verified from the defini- 
tions in Section 2. The integral was numerically evalu- 50 
ated for different estimators and delays. FIG. 6 shows 
the loop bandwidth as a function of the estimator time 
constant and N for type II and type III (second order 
and third order) estimators. Observe that the introduc- 
tion of extra delays in the loop increases the noise band- 55 
width, but the effect of the delays is less significant as 
the estimator time constant increases. 

(3) Stability and Gain Margin 

The closed loop root-loci for the designed DPLL 60 
will now be discussed, recalling the definition of the 
open loop transfer function in equation (12). The open 
loop transfer function G(z) has zeros at the zeros of 
H(z) and poles at the zeros of Ej(z), thus, with increase 
of the phase detector gain A, the poles of P(z) move 65 
from the zeros of Ei(z) to the zeros of H(z). 

FIGS, la, lb , 1c and Id show root-locus plots for the 
type II and type III loops, for the case N=2, The trans- 


it) 

fer function Ei(z) has three zeros, two at z— 1 and one 
at z— — 2(1 — ct). Observe in FIG. la that at low gain, 
the poles of the closed loop transfer function start at 
these points; then, as the gain increases, they move 
within the unit circle. Eventually, with higher gains, the 
poles of P(z) go outside the unit circle and stability is 
lost. The root locus for the third order estimator illus- 
trated in FIG. lb has a similar form, except that at very 
low gain, roots occur outside the unit circle, near z = 1. 

As shown previously, the designed approach can be 
marginally stable due to cancellation of a pole of S(z) by 
a zero of the NCO. To assure stability, this pole at 
z= — 1 should be slightly shifted to the inside of the unit 
circle. Effect of this shift in pole location is demon- 
strated in FIG. 7c and FIG. Id . These figures repeat the 
cases given in FIG. la and FIG. lb, with new pole 
location at z= — 0.98. There is an extra root locus 
branch from z= — 1.0 toward Z=— 0.98, with negligi- 
ble effect on the rest of the root loci. In fact, since the 
whole root locus branch is inside the unit circle, an 
arbitrarily small shift in pole location is sufficient to 
guarantee loop stability. The small shift in pole location 
has a negligible effect on the other parameters of inter- 
est: dynamic phase error, gain margin, and noise band- 
width. 

Locations of the poles of the closed loop transfer 
function at the operating point are of special impor- 
tance. These poles are commonly selected to be on the 
real axis so to avoid instability and oscillations when the 
gain A changes slightly. At the operating point, A=l, 
the closed loop transfer function is H(z). For the estima- 
tor-predictor in our example illustrated, H(z) has two or 
three poles at z=0 and two or three co-located poles at 
z—a for the second and third order estimator, respec- 
tively. This feature, real closed loop transfer function 
poles at the operating point, satisfies the stability criteria 
often cited by PLL designers. 

Gain margin measures the effect of the phase detector 
gain, A, on the stability of the closed loop. It is defined 
as the ratio of the critical gain, i.e., gain that forces the 
poles of the transfer function outside the unit circle, to 
the nominal unity gain. FIG. 8 shows the relationship 
between upper gain margin and normalized bandwidth. 
Here, as the loop becomes wider, the gain margin de- 
creases. Similar results are presented in FIG. 9 for the 
lower gain margin of the type III loop. Note that the 
lower gain margin for the type II loop is infinity. These 
figures are useful for quick design evaluation e.g., find- 
ing a favorable trade-off between any two variables. 

Similar results can be obtained for a tradeoff bet ween 
dyamic errors and B^T. Since dynamic phase error is 
approximately inversely proportional to the second or 
third power of the normalized bandwidth, we define 
phase error coefficients, C<j>2 and C<p3, so that: 


r c * 2 i 

Type II: £00 — I acceleration 

r Q* "| 3 

Type III: Eco « [ — J jerk 
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These coefficients exhibit lower variation with BlT, 
than the dynamic phase error. FIG. 10 shows the phase 
coefficients as function of normalized bandwidths. The 
phase coefficients, for narrow bandwidths, are less than 
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1 for the type II loop and between 1 and 2 for the type 

III loop, and they increase slowly with increased band- H(z)-C(z)D(z)z- N 

widths. 


Using an estimation theory approach, design of this 
new DPLL can be accomplished in a systematic way 5 
rather than by trial and error. The examples presented 
here are useful for many applications, and the resulting 
design curves may thus be directly applied. Still other 
examples and variations will occur to those skilled in 
the art. Consequently, the claims appended hereto are 10 
intended to encompass other equivalent examples and 
variations. 

What is claimed is: 

1. A digital phase-lock loop comprising a linear esti- 
mator means for recursive estimates of phase, fre- i5 
quency, and higher order derivatives, a predictor means 
for compensating transport lag inherent in said loop, a 
phase detector means with gain for comparing an input 
signal Q(z) with the output Q(z) of said predictor with 
a delay inherent in said loop, and a means for adding the 2 
output of said predictor having said inherent delay z _Ar 
to the output <j>(z) of said detector, said loop having a 
closed-loop transfer function of the form 


where C(z) is an estimator transfer function, D(z) is a 
predictor function and z-^is said transport lag in units 
of loop update time inherent in said loop. 

2. A digital phase-lock loop as defined in claim 1 
wherein said estimator means and said predictor means 
with inherent transport delay are comprised of a digital 
filter having a transfer function S(z) and a numerically 
controlled oscillator having a transfer function Q(z), 
where said transfer function S(z) satisfies: 

S( 2 ) — . £( z ) D( z ) 2 -n 

() m z ■ 

3. A digital phase-lock loop as defined in claim 2 
wherein said numerically controlled oscillator, said 
digital filter and said adding means are implemented by 
a programmed digital computer, and a separate numeri- 
cal controlled oscillator having a transfer function Q(z) 
is provided between the output of said digital filter and 
said phase detection means. 

* * * * * 
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